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We extend the recently developed non-gaussian ther- 
modynamic formalism Q of a (presumably strongly tur- 
bulent) non-Markovian medium to its most general form 
that allows for the formulation of a consistent thermo- 
dynamic theory. All thermodynamic functions, includ- 
ing the definition of the temperature, are shown to be 
meaningful. The thermodynamic potential from which 
all relevant physical information in equilibrium can be 
extracted, is defined consistently. The most important 
findings are the following two: (1) The temperature is 
defined exactly in the same way as in classical statistical 
mechanics as the derivative of the energy with respect to 
the entropy at constant volume. (2) Observables are de- 
fined in the same way as in Boltzmannian statistics as the 
linear averages of the new equilibrium distribution func- 
tion. This lets us conclude that the new state is a real 
thermodynamic equilibrium in systems capable of strong 
turbulence with the new distribution function replacing 
the Boltzmann distribution in such systems. We discuss 
the ideal gas, find the equation of state, and derive the 
specific heat and adiabatic exponent for such a gas. We 
also derive the new Gibbsian distribution of states. Fi- 
nally we discuss the physical reasons for the development 
of such states and the observable properties of the new 
distribution function. 

05.20.-y, 05.70.Ce, 51.10.+y, 52.25.Dg, 52.35.Ra, 52.65.Ff, 
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I. INTRODUCTION 

In a recent paper we developed the kinetic theory of 
a coUisionless (presumably non-Markovian) equilibrium 
state of a system of N particles undergoing strongly tur- 
bulent interactions. Our motivation was to investigate 
what kind of statistical mechanics described an equilib- 
rium state of fully developed stationary turbulence, if 
it existed at all. Evidence of the possibility of such a 
description comes from observations of Levy flights and 
started in the twentieth of this century when Richardson 
found his turbulent spectral law that was fundamentally 
different from Kolmogorov's law of ordinary spectral be- 
haviour in turbulence (for a review of the history see, 
e.g., |). 

There had been earlier attempts to construct other 
than Maxwell-Boltzmann probability distributions al- 



ready in the past century and continued until today, 
mostly as attempts to describe the thermodynamics of 
extreme situations ||. These attempts culminated in a 
purely mathematical extension of Boltzmann's defini- 
tion of entropy by Renyi ^,|^. These extensions, though 
widely used in chaotic dynamics in order to infer about 
multi-fractal behaviour, have not been given any physi- 
cal justification yet. A certain mathematical simplifica- 
tion of Renyi's original proposal has subsequently been 
suggested || that seemed to have applicability to some 
range of physical problems. In particular, some kind of 
thermodynamics was constructed (for the most lucid pre- 
sentation see, e.g., @). However, this theory refers to an 
unusual prescription of constructing physical observables 
that is not in agreement with conventional physics. 

Our basic assumption was that the stationary turbu- 
lence we were going to describe was not describable by 
weak turbulence theory. This assumption implied that 
any equilibrium state the system might have achieved 
could not be constructed by the means of perturbation 
technique, i.e. no small expansion parameter exists. Sys- 
tems of this kind are subject to non-perturbation tech- 
niques. In some of those systems the transition from the 
original state towards turbulence proceeds through crit- 
icality. Such systems are critical and must be treated 
by renormalisation group methods (cf., e.g., [|||]). We 
did not explicate about this point but instead asked for 
the properties of the corresponding equilibria, propos- 
ing that the system had actually settled into a critical 
and turbulent intermediate equilibrium. We found that 
these equilibrium states were describable by a new non- 
Boltzmannian collision integral 



dr-^dOG T [12], 



(1) 



where ctt is the (turbulent) collisional cross-section, and 

G T [12] = g[f(l')]g[f(2')} - <?[/(l)M/(2)] (2) 

is the correlation functional of the one-particle distribu- 
tion functions / after (primed) and before (unprimed) 
the interaction. Gt[12] itself was found to be the prod- 
uct of functionals g[f] each depending on the distribution 
function of one family of particles only. (The index T in- 
dicates that the systems are in a state of about stationary 
turbulence.) The turbulent state in that they are found 
is reached on passing through critical points after having 
entered a highly nonlinear phase. The transition to the 
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turbulence is not known in detail, but it can be assumed 
that during the critical phase the systems evolves into 
all scales such that the scales cannot be separated any- 
more. It becomes essentially scale-invariant, and any per- 
turbation theory breaks down when all scales are highly 
correlated. In pj it had been demonstrated that the func- 
tionals g[f] actually consist of infinite products of corre- 
lations suggesting that this interpretation is close to the 
truth. 

The breakdown of perturbation theory in the scale- 
invariant state suggests that a microscopic approach to 
the problem will have to refer to renormalization group 
techniques. We have shown in jjj, however, that it is 
possible to describe the stationary equilibrium state of 
the system without the need to develop the microscopic 
phase-transition theory in detail by the methods of sta- 
tistical mechanics. This has been done in close analogy to 
Boltzmann's kinetic theoretical approach. Detailed bal- 
ance then requires that 



lng[f] = -/3(ep - m) 



(3) 



where e p = p 2 /2m is the energy of a particle of momen- 
tum p and mass m, and /3, y, are two arbitrary constants 
(playing the role of Lagrangean multipliers) which have 
been suggested jij to correspond to the kinetic temper- 
ature and chemical potential of the system, respectively. 
Below this suggestion will be proved in full strength. 

If the system is described by Equation (|l|) and is in 
the assumed highly nonlinear multi-scale turbulent equi- 
librium, then one can introduce ^] a control parameter 
k such that 



lim g[f(n,e p )] -> / B (e P ) 



(4) 



reproduces the Boltzmann distribution function /b, and 
Gt becomes the ordinary Boltzmann collision functional. 
In |Q a particular functional g[f) was found and it was 
demonstrated that this functional actually described a 
thermal equilibrium state nicely satisfying an iJ-theorem 
and permitting for a new expression for the turbulent en- 
tropy St. One should note that this new mathematical 
expression for the entropy does not attach any new phys- 
ical interpretation to the entropy. As in the Boltzmann 
case, entropy describes the amount of irreversible disor- 
der in the system. The new expression for the entropy 
merely means that in the turbulent scale-invariant state 
of the system the increase of disorder is calculated in a 
different way than in the conventional Boltzmann state. 



where {k,£} € R, and £ is an arbitrary fixed real num- 
ber (R is the space of real numbers). The advantage of 
introducing £ will become clear below. It is then easy to 
demonstrate that the condition (||) is satisfied for any 
arbitrary fixed I ^ oo. 

The equilibrium distribution function /(k, e p ) can be 
constructed using the two equations and (0). The 
most general distribution function is found to be 



/3/z /3ep 

K K 



(6) 



Clearly, it is a function of particle momentum p through 
the particle energy e p = p 2 /2m and of the two param- 
eters, k and I, respectively. It is then easy to show in 
parallel to jjj that fe satisfies the following new gener- 
alised turbulent entropy relation referred to in the above 
discussion: 



fp 
/i 3 



fi(e p ,K)lng[fi(e p ,K)}. (7) 



This entropy is concave and moreover is super-additive, 
meaning that the entropy of two independent systems 
is larger than the sum of the individual entropies of the 
two systems, an interesting property of such states. (Note 
that, in accord with physical intuition and requirement, 
St,£ can never become sub-additive. Entropy in a single 
closed system or in a collection of closed systems will al- 
ways grow as disorder cannot be extinguished by adding 
other disorder.) That this is true can be shown along 
the same lines as in |l[ , for the introduction of the fixed 
number I does not introduce any change in the analysis. 
Hence, ft is an actual thermodynamic equilibrium distri- 
bution that replaces the Boltzmann distribution /b under 
K-conditions. All the discussion of 0j can be applied to 
it. 

In the following we develop the corresponding thermo- 
dynamics and show that it requires the choice I = 1. 
Therefore, the correct thermodynamic equilibrium one- 
particle distribution function of the K-gas is given by 



f(e p ,n) =1 + — £ 

K K 



(8) 



where, for convenience, we suppressed the index I = 1. 



III. THERMODYNAMIC RELATIONS 



II. GENERALISATION 

We now generalise the functional g[f] to its most general 
form 

g[f( K , e p )] = exp {k [l - r 1/{K+e \n, e p )] } (5) 



In this section we define the basic thermodynamic func- 
tions and demonstrate that a K-gas is a system that 
though behaving in a special way is nevertheless in com- 
plete accord with the fundamental thermodynamic rela- 
tions. 
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A. Thermodynamic Potentials 

All macroscopic thermodynamic information about a sys- 
tem in contact with the external world is contained in the 
thermodynamic potential Q(/3, V, /i) which is a function 
of the temperature variable /?, the volume of the system, 
V, and the chemical potential \x (or particle number N). 
We define Qt{P, V, /i) of the K-gas by 



Qt(/3,V» 



d 3 p 



Pv . /?e P 



(9) 



Since we assume that thermodynamics should provide a 
valid description of the macrostate of the K-gas in order 
to be in accord with conventional physics, the relation 
between the thermodynamic potential and other ther- 
modynamic functions is given by 



Q T = F T - fiN. 



(10) 



Here Ft is the Helmholtz free energy of the K-gas, and 
Qt is the difference between F and the product of the av- 
erage particle number and the chemical potential. Hence 
the free energy is given by 



F T = Qt + fJ.N = - 



V 



d 3 p 



Pa 

K 



■fiN. 
(11) 



We now show that this definition is consistent with the 
Boltzmann limit. Indeed, taking the limit k — > oo in the 
integral we immediately find that 



F = ^ J IF ° XP h/?(£p " ^ )] + ^ 



(12) 



which coincides with the Boltzmann-Helmholtz free en- 
ergy. Moreover, this expression identifies fi with the 
chemical potential, and [3 = with the inverse ki- 

netic temperature. The latter expression will be proved 
explicitly and in full generality for arbitrary k below. 



B. Average Energy and Number Density 

In conventional statistical mechanics, the total mean en- 
ergy E is calculated from the Helmholtz free energy in 
the following way: 



E 



d 



dp 
d_ 

~d/3 



(PF) 



yj^-exp [-P{e p - /i)] - fiN 



(13) 



In the new thermodynamics we replace F with Ft to 
obtain 



d_ 

~d/3 



V 



d 3 p 



(14) 



Carrying out the partial differentiation, this yields the 
following expression for the energy 



E = V 



V 



d 3 p 



h 3 (1 - pn/n + j3e p /K) K 



+i 



d 3 p 



e p f(e p ,K) 



(15) 



Note that the second part of this equation is just the 
correct physical definition of the average energy of the 
system as the integral over phase space of the distribution 
function, Equation (||), as is required by the commonly 
used definition of an observable in statistical mechanics 
and kinetic theory. In this respect, our theory is thus 
consistent with common statistical physics and kinetic 
theory. It does not require any different kind of averaging 
in order to calculate the observables as the physically 
relevant quantities. This implies that this theory is also 
in accord with the fundamental kinetic BBGKY theory 
(cf., e.g., @, chp. 2). 

In order to be consistent with the previous sections, 
Equation ([l5] suggest that we must identify I = 1, which 
justifies our previous choice for i. It would of course be 
possible to choose any arbitrary I in the definition of Q 
and to adjust the exponent of / appropriately, but such 
an action would introduce some unnecessary arbitrariness 
that would result in a simple re-scaling of n. 

The average particle number N is then given as the ze- 
roth order moment of the distribution function /(e p ,«;). 
This can be shown to be the negative of the partial 
derivative of the above thermodynamic potential Qt 



^=v7^[i-^+^ 

J h 6 K K 

[OQt 



-(k+1) 



(16) 



The average density of the K-gas is correspondingly ob- 
tained as 



N 
V 



d^p 

h 3 



(17) 



On the other hand, it must be required that the chemical 
potential of the K-gas is the partial derivative of the free 
energy Ft with respect to the average particle number 



/i = {dF T /dN) 0v . 



(18) 



This equation is conventionally used to express fi through 
N. 
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C. Entropy 

In order to find the entropy relation, we form the 
following derivative of the thermodynamic potential 
kB(dQ,T/df3~ 1 )v fJ ,- It can be shown that this derivative 
can be written as 



3(1//?) 



= -kB0[E -fiN- Q T ] 



E). 



(19) 



The last expression is just the negative of the entropy St 
and coincides with the following definition 

t d 3 p 

St = -k B V / — 3- f{e p ,K)bxg[f(e p ,K)] 

= -k B V J ^/(c P ,«)[l-r 1/(,,+1) (cp,/c)] (20) 

when replacing the distribution function inside the inte- 
gral with Equation (||) and using the above derived rep- 
resentations for the average energy ( |l5| ) , particle number 
(|l6|), and the thermodynamic potential (0). 

Clearly this entropy expression has a structure differ- 
ent from the ordinary Boltzmann definition. It contains 
the logarithm of the correlation functional g[f] in place 
of the logarithm of the distribution function itself. This 
fact implies that in the K-state of the gas it is the cor- 
relations that contribute most to the entropy. While in 
the final state of the system when the interactions be- 
come purely stochastic and the gas settles into its ther- 
mal death, g — > /b, and the entropy assumes its classical 
representation. 



D. Consequence I: Definition of Temperature 

As an important application of the above definition of the 
entropy St of the K-gas we now derive the expression for 
the thermodynamic temperature. From classical thermo- 
dynamics it is known that the only consistent definition 
of the temperature is given in the form of a derivative of 
the entropy: 



1 

T 



dS_\ 



(21) 



When we use the expression (||) for the distribution func- 
tion in the definition of the entropy (ECh we recover that 



S T = k B p(E- tiN), 



(22) 



which is nothing else thatn a rearranged version of (19). 
Then taking the partial derivative with respect to the av- 
erage energy E we immediately identify the temperature 
T as in conventional thermodynamics with the inverse of 
the Lagrangean multiplier (3 



J3 = feT - 



(23) 



This very important relation proves that the tempera- 
ture T of the K-gas is defined exactly in the same manner 
as the temperature of the classical Boltzmann gas. In this 
way it renders all other definitions of the temperature of 
the K-gas used in the literature invalid. Those definitions 
still contained dependencies on k (cf., e.g., |l3|-|l6[] and 
elsewhere) . These K-dependencies turn out to be unphys- 
ical. They are the result of a naiv use of the so-called ex- 
perimentally determined '^-distributions' in calculating 
a formal expression for the temperature. The physically 
correct distribution function Equation (H) resembles the 
/^distribution, but it contains the non-vanishing chemi- 
cal potential. It is only this complete distribution func- 
tion that leads to a correct thermodynamics, with the 
temperature T being defined as a physical quantity by 
its thermodynamic definition, Equations (||l]) and (|23|). 
The temperature in this sense is a measure of the state 
of the system. It is a parameter that characterises the 
gas. It is not a measure of the mean energy. Only in the 
Boltzmann gas the two, mean energy and temperature, 
measured in energy units, are related in a simple way. 
Below we are going to derive the relation valid in the 
K-gas. 

The above derivation of the temperature thus justifies 
the use of the constant /3 as the unambiguous and only 
measure of the thermodynamic temperature of any (tur- 
bulent) K-system and shows that in such a system the 
above definition of the entropy consistently replaces the 
Boltzmann definition. 



E. Consequence II: Equation of State 

In order to complete the set of thermodynamic relations 
we may now construct the equation of state taking the 
partial derivative of Qt with respect to the volume V . 
This procedure yields the pressure 



P 



dQr 
dV 

d 3 p 



P» 



/3(i (3e p 



(24) 



Rewriting this expression with the help of Equation (|^) 
we find the following important relation 



PV = -Q T (/3,V,n). 



(25) 



This is the fundamental equation of state of a (turbulent) 
K-gas. It shows that such gases possess complicated equa- 
tions of state. Such behaviour has been expected from the 
very beginning, because the presence of the long-range 
correlations should become manifest in the average prop- 
erties of the gas as well. Even an ideal K-gas turns out to 
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have a non-simple equation of state, as will be demon- 
strated below. 

The remaining first and second order thermodynamic 
relations can all be obtained from the previous relations 
and will not be given here. 



IV. REVIEW OF IDEAL k-GAS PROPERTIES 

In our previous paper B we derived some of the relations 
for an ideal gas. Here, because of the precise definition of 
the thermodynamic potential and the above consistency 
proof of the thermodynamic relations we can consider- 
ably simplify the expressions given before. Moreover, the 
formulas given below correct them for thermodynamic 
consistency. 



A. Chemical Potential 

We start with the average particle number N of an ideal 
re-gas. This number is given as the phase-space integral 
over the ideal gas distribution function /(e p , re) Equation 



( 



N = V 



T(k- 1/2) 
T(k+ 1) 

o \ 3/2 



h 2 (3 



fa 

K 



"(k-1/2) 



(26) 



It is assumed throughout thermodynamics that the av- 
erage number density n = N/V would be known. Hence, 
Equation (^6|) is basically an equation for the chemical 
potential /i of the ideal re-gas. Inverting (^) and intro- 
ducing the 'quantum re-density' 



/ 2irmK, \ 3 ^ 2 

nqK = \ W ) 

we obtain for the chemical potential 



fa 



n qK T(k - 1/2) 



r(« + i) 



l/(«-l/2) 



(27) 



(28) 



The right-hand side of this expression is always positive, 
and hence either the chemical potential is negative /i < 
or fa/n < 1. It is not too difficult to demonstrate 
that in the limit re — > oo this expression reproduces the 
classical chemical potential of an ideal gas |10| with n q — 
(2nm/h 2 /3) 3 / 2 the 'quantum density' of the classical gas. 



B. Mean Energy 

Also, carrying out the integration in Equation (|l5|), the 
average energy of the ideal re-gas follows as 



E 



3N 



2/3 « -3/2 

3iV k 
2/3 re - 3/2 



fa 



n qK T( K - 1/2) 
n T(k + 1) 



l/(«-l/2) 



(29) 



With the above identification of j3 = 1/k^T one imme- 
diately realises that the limit re — > oo reproduces the 
classical well-known result (cf., e.g., p0| , p. 77) 



E 
N 



-knT 



(30) 



for the energy per particle. As is obvious from Equation 
( p9| ) , in the re-gas this last relation between the mean en- 
ergy and the temperature becomes much more involved. 
Moreover, it is also obvious that the above definition of 
the energy requires that 



3 

2' 



(31) 



a condition that is consistent with a similar one derived 
in our previous publication JjJ. 

C. Ideal Gas Equation of State 

We now calculate the pressure of the ideal re-gas in 
order to determine its equation of state. From Equation 
(p5|) we obtain 



PV 



N 



re- 
N 



3/2 



1 



fa 
re 



(3 re - 3/2 



n qK T(re - 1/2) 
n T(re + 1) 



l/(«-l/2) 



(32) 



This equation of state becomes the ordinary ideal gas 
equation of state PV — Nk-gT only in the Boltzmann 
limit re — > oo. The ideal re-gas behaves differently, pos- 
sessing a much more complicate equation of state. This 
seems reasonable because, as mentioned above, the fact 
that the re-gas is a highly correlated many-body system, 
which we assume is in an evolved turbulent state, must 
in the first place become obvious in its equation of state. 
Such a gas should exhibit a behaviour that differens from 
that of an ordinary laminar ideal gas. Nevertheless, how- 
ever, is is then both surprising and satisfactory that, com- 
paring Equations (^) and (|32|), one recovers the basic 
thermodynamic relation between the pressure and aver- 
age energy, 



PV = -E, 



(33) 



as is valid also in ordinary ideal gases. The unbroken va- 
lidity of this relation even in the re-gas is interesting and 
is an important finding. Its physical meaning is that it 
relates the mean energy in a simple geometrical way to 
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the pressure of the K-gas. This retains the physical in- 
terpretation that the pressure is the geometrical effect of 
the average motion of the constituents of the gas. The 
temperature assumes quite a different meaning then as 
an independent parameter measuring the irreversibility 
of the system in place of its thermal mean energy. The 
average energy of the K-gas is thus not a simple measure 
of the temperature T of the K-gas. The average energy 
per particle E/N in a K-gas is not anymore a simple frac- 
tion of the kinetic temperature. Only in the Boltzmann 
gas become both, mean energy and temperature, simple 
equivalents. 



D. Specific Heat and Adiabatic Index 

In order to complete the discussion of the ideal gas, we 
calculate the specific heat Cy at constant volume V from 



CV — -/3 2 fcB ( TT77 
Since E is known, it is simple matter to find 



(34) 



C v = ik B N 



3/2 
T(k- 1/2) 



nl/(K-l/2) 



T(k + 1) 



(35) 



Again, in the limit k —> oo, this expression converges to 
the ideal gas value (3/2)/cb^V as the quantity in the curly 
brackets becomes just 



lim (2TTm/h 2 l3) 3/i2K ^ 1) 



1. 



(36) 



However, as expected, Cy in the K-gas is not a constant 
(or otherwise an extensive function of the particle num- 
ber). It depends on both, density and temperature. This 
complication suggests that the adiabatic exponent 7 K of 
the /t-gas will as well not be a simple constant, as is the 
case for the ideal Boltzmann gas. 

Calculation of the specific heat Cp at constant pressure 
turns out to be more involved. We can, however, take 
advantage of the general thermodynamic formula 



Cp — Cy — T 



(d(3/dV)i 



(37) 



relating Cp and Cy (cf., e.g., jlXJ). This is possible be- 
cause we have already demonstrated that the thermody- 
namics of the K-gas is formally identical with that of the 
ordinary Boltzmann gas. 



Making extensive use of the equations of state (32), 
energy (p9|), and specific heat at constant volume ( |35| ) 
we obtain for the difference (|37|) of the two specific heats 



Cp — C\ 



knN- 



k- 3/2 

TV r 1\ -|l/(«-l/2) 



n qK T(K- 1/2) 



(38) 



Inspection of this equation suggests that in a dilute k- 
gas where the density n = N/V <C n qK is much less 
than the quantum density, the difference between the two 
specific heats will become small. On the other hand, close 
to maximum correlation (at k — ► 3/2) this difference can 
become large. These are two further distinctions between 
ordinary and K-gases, respectively. 

The most interesting quantity is, however, the index 
of adiabaticity, 7 = Cp/Cy. Dividing by Cy in (|3q ) one 
finds that in the K-gas this ratio becomes 



7« 



1 

1+ 3 



t(k + i) 



n,„r(K-l/2) 



1/(k-1/2) 



n qK T( K - 1/2) 



T(k+ 1) 



l/(«-l/2) 



Indeed, for k 



one can show that 



lim 7 K 



3' 



(39) 



(40) 



which is in accord with the thermodynamics of an ordi- 
nary (three-dimensional) gas. But, for k < 00, the adia- 
batic index may deviate considerably from this value. In 
particular, when the term in the brackets containing the 
quantum density is sufficiently large compared to 1/2, 



7« 



'♦5 



r( K + 1) 



n qK r( K - 1/2) 



2/(k-1/2) 



(41) 



This value is clearly smaller than 5/3 and in the limit of 
a very dilute gas approaches unity. In particular, taking 



3/2, 



73/2 ~ 1 + 



18 \n 



(42) 



This is very close to unity for small ratios n/n q . For clas- 
sical K-gases we can therefore conclude that the adiabatic 
index falls into the intervall 



K Ik < 7 = g 



(43) 



with the right-hand side holding for the ordinary gas. 
On the other hand, there is the possibility that in very 



dense K-gases with n 



' QK 



the second term in Equation 



(|39j) becomes negative. In such a situation, values 7 K < 1 
would even arise. For instance, assuming n — * n q , one 
finds that 



lim 



73/2 



V2tt 
3%/3-2V2 



< 



(44) 



G 



becomes negative. Of course, this extreme case is phys- 
ically impossible, and is a result of using the minimum 
value for k, and n — n q at the same time. This discussion 
shows, however, that for sufficiently high densities and 
small k the adiabatic index becomes small. Investigation 
of these cases requires the development of the quantum 
theoretical extension of the present theory. This will be 
given elsewhere p|. 

The above result on the reduction of the adiabatic in- 
dex is exciting. Reduction of 7 physically implies that the 
number of degrees of freedom increases. Remember that, 
classically, 7 is a measure of the internal freedom of the 
system and is microscopically related to the degrees d of 
freedom by 



7 = (d + 2)/d. 



(45) 



Equation ( f43| ) then suggests that for K min < n < 00, the 
number d increases by a large amount. For d > 2, the 
index 7 quickly approaches the value 7 — > 1. In the k- 
gas, 7 K may not necessarily be related to d in the same 
simple way (45) as in the classical gas. It is possible that 
the relation between 7 K and d will be more involved. In 
order to obtain the actual dependence, one must develop 
the microscopic theory of the interactions in a multi-scale 
medium which is outside our present reaches. Neverthe- 
less, the decrease in 7 provides a strong argument for the 
validity of the initial assumption underlying our theory. 
This assumption was that the K-gas is in a multi-scale 
state with an enormously large number of degrees of free- 
dom, and that these scales are all correlated in a way as is 
expected for fully developed turbulence. Media like this 
are believed to be scale-invariant and should microscop- 
ically be treated by renormalisation-group methods |J . 



in the square brackets on the right-hand side of Equa- 
tion ( [46| ) is large, we still arrive at the surprising result 
that the condition Equation ( f47j ) is still exactly satisfied. 
Hence, only the intermediate regime when 



" , qK 

n 



r(« + i) 



!/2r( K - 1/2) 



(48) 



the adiabatic relation between volume and temperature 
changes. This practical independence of the condition of 
isentropy on the value of kappa is a purely geometrical 
effect that tells that a sufficiently fast expansion of any 
medium should cause cooling. It is only reasonable that 
this behaviour does not halt in front of turbulent or scale- 
invariant systems. With the help of (^7|) the two remain- 
ing isentropic relations are obtained from the equation of 
state as 



pT -b/2 = congt; 
PV 5/3 = const. 



(49) 
(50) 



It is remarkable that these conditions exactly coincide 
with those for the isentropic processes in a classical Boltz- 
mann gas. We thus learn from this agreement that isen- 
tropy results in simple geometric behaviour of both Boltz- 
mann and K-gases. Since the thermodynamic differential 
equations remain valid for the K-gas, they can in prin- 
ciple be used to express these exponents through j K . In 
a K-gas, however, the ratio of specific heat has no direct 
relation to the isentropic processes as is the case in an 
ordinary ideal gas. 



F. Volume Coefficients 



E. Isentropic Process 

With the adiabatic index at hand we can investigate the 
isentropic (adiabatic) evolution of the K-gas. To this end 
we need the explicit expression for the ideal gas entropy 



St = kftV 



2irKm\ 3/2 r(K-3/2) 




h 2 P J 
n qK T(k - 



1/2) 



T(k+1) 



1/(k-1/2)- 



(46) 



Isentropic processes leave the entropy constant. To first 
approximation this requires that 



VT 3/2 = const, 



(47) 



a condition that is identical to the condition for adia- 
baticity in the thermodynamics of ordinary gases (cf., 
e.g., 0|). For K-gases, the second term in the curly brack- 
ets in Equation ( fl6| ) seems to introduce a further com- 
plication. However, in the particular case when the term 



In order to conclude this section we finally derive the 
coefficients of thermal expansion and compressibility, re- 
spectively. The first is given by the well-know formula 



v\df) P , 



(51) 



It can be easily calculated from the equation of state 
( p2] ) of the K-gas. Replacing the chemical potential in 
the final expression after having performed the relevant 
differentiations, we find that 



1 
T 



r(/s + i) 



i qK r(«- 1/2) 



1/(k-1/2) 



(52) 



retains its inverse proportionality to the temperature T 
while otherwise being small for dilute media. The thermal 
compressibility 



V \dP) 



(53) 
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can most easily be obtained from Equation 
ing the relation 



C P -C V = TV- 



*cxp 

Kt"' 



exploit- 



(54) 



This relation together with Equation ( |52| ) leads to the 
following useful result 



V re -3/2 
iVfc B T re 

?! r(« + i) 



n ?K r(« - 1/2) 



l/(«-l/2) 



(55) 



This completes our account of the stationary properties 
of the ideal re-gas. 

V. FLUCTUATIONS 

Because it is a system in thermal equilibrium the re-gas 
is capable of thermal fluctuations as well. Such fluctua- 
tions are known to be the root-mean-square amplitudes 
of the oscillations of the various macroscopic quantities 
around the thermal equilibrium state. Independent of the 
very particular form of the equilibrium distribution func- 
tion the fluctuation amplitudes average out when aver- 
aged over sufficiently long times or spatial scales. Under 
re-conditions this restriction poses caution on the defi- 
nition of fluctuations. Times can be no longer than the 
inverse binary collision time as it is our philosophy that 
the collisionless turbulent scale-invariant state refers to 
times shorter than the binary collision time 1/V C (cf., [EJ). 
Hence, very slow oscillations may not average out even 
linearly and may survive the entire scale- invariant regime 
until they enter the binary collisional state when they be- 
come ultimately depleted. Moreover, the assumption of 
scale-invariance (or self-similarity) also implies that very 
long oscillations may survive. With these restrictions in 
mind the rms fluctuation amplitude of a quantity A is 
defined in the ordinary way as 



(AA) 2 ms = (A 2 ) (A)' 



(56) 



where the average is understood as the ensemble aver- 
age or expectation value. Our discussion presented in the 
previous sections has shown that calculating expectation 
values by linear averaging over the distribution function 
is the appropriate way to stay in accord with the require- 
ment of being able to define a consistent thermodynamic 
theory in the re-regime of the evolution of the system. 
Here we are interested in providing the expressions for 
the most fundamental rms fluctuation amplitudes. 



A. Particle Number Fluctuations 

The particle number fluctuation cannot be calculated di- 
rectly from the zeroth order moment of the distribution 



function. Remembering, however, that the particle num- 
ber is closely related to the chemical potential the general 
thermodynamic expression for the number fluctuation is 
(cf., e.g., M p. 152) 



(AN) 2 ms = k B TV 



d 2 P 



(57) 



The second-order derivative of the pressure on the right- 
hand side of this equation can be shown to be equiv- 
alent to 



8 2 P 
9/i 2 



TV 3 
V 3 



dP 



_d(V/N) 



(58) 



Using the ideal gas equation ( |32| ) for the pressure, we 
obtain for the rms particle number fluctuation 



(AN) 



2 

rms 



3/2 



N 2 



Nk 

n r(re - 



n, K r(re-l/2) 



1/0-1/2) 



(59) 



The rms number fluctuation increases with N while the 
relative fluctuation in particle number decreases as N in- 
creases. Note that if we would have used the average ex- 
pressions to calculate the number fluctuation, we would 
have found a zero fluctuation in number. Hence the effect 
is small in a dilute g it should be. Moreover, com- 
paring the fluctuation with the compressibility (^5|) one 
observes that the number fluctuation is proportional to 
Ktk as has been expected since this is a particular case 
of the fluctuation-dissipation theorem holding for small 
fluctuations (cf., e.g., chpt. 8). 



B. Energy Fluctuations 

It is more interesting to determine the fluctuation of 
the mean energy. It is most convenient to simply calcu- 
late the two expectation values in Equation (|56l ) with 
A = E and to take their difference. It is, however, much 
simpler to use the general thermodynamic expression for 
the energy fluctuation 



(AE)l 



-(dE/d/3) 



(60) 



in order to obtain the average fluctuation of the internal 
mean energy 



(A£)? ms = 



2 \k- 3/2 



N(k B Tf 



(61) 



The rms fluctuation amplitude is proportional to the 
thermal energy k B T and increases as the root of the 
particle number. Its absolute value also increases with 
re approaching re m i n . The smaller re the higher is the ab- 
solute level of the fluctuations in energy. On the other 
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hand, when dividing the rms amplitude by the mean en- 
ergy Equation (^9|) one finds that at the same time the 
relative energy fluctuation 



E 2 



1 k — 3/2 
N K 



(62) 



decreases because the mean energy increases faster than 
the fluctuation. This can be easily understood as the in- 
creasing effect of storage of energy in the high-energy 
tail of the distribution near n mul . The relative amount 
of energy stored in the fluctuations produced by the tail 
becomes less than the energy itself. Nevertheless, though 
the energy diverges faster than the fluctuation energy, 
the gas contains fluctuations of large amplitude that is 
an expression of its turbulent nature and the appearance 
of many scales. 



VI. PROPERTIES OF THE EQUILIBRIUM 
DISTRIBUTION 



In this section we return to the equilibrium distribu- 
tion function /(c p ,k) and investigate its most obvious 
properties. Because distribution functions nowadays have 
become measurable quantities (though not observables 
in the strict physical sense) their shape and behaviour 
provide a direct link to the physical properties of the 
medium. Examples of so-called ^-distributions observed 
in space plasmas have been given continuously for about 
thirty years (cf., |jl9|-p3"|), with the most extended ob- 
servations being obtained in the Earth's magnetospheric 
tail |2l|]22| |. Similar distributions are believed to be re- 
lated to Levy flight dynamics [|| and have been used in 
an attempt to interpret cross- field diffusion pij . But the 
physics of such K-distributions has not been understood 
for a long time. Attempts of their explanation trace back 
to various types of solutions of the Fokker-Planck equa- 
tion for particular interactions |p5 - p7[ . The general ther- 
modynamic theory developed in |1 1 and in the present pa- 
per is the first consistent physical theory of these families 
of distribution functions, aside from another approximate 
attempt p8[ that was based on so-called non-extensive 
thermodynamics 

In the following we discuss the information that can 
be extracted from observation of generalised-Lorentzian 
distribution functions. 



A. Velocity Distribution 

The distribution function /(e p ,K) Equation^ has been 
given in terms of the particle momentum p. It belongs to 
the family of generalized-Lorentzian functions. Distribu- 
tion functions like this one possess the property that not 
all infinitely many moments of the distribution can exist. 



For given n, the moments of order r > 2k — 1 start di- 
verging. This fact restricts the value of k to those values 
that satisfy the condition 



k > (r + l)/2. 



(63) 



For the mean energy to be a physically measurable quan- 
tity this expression sets the limit for k > 3/2. At the 
current state it is difficult to discuss what will happen to 
the higher moments and in which way the theory can be 
'renormalised' if at all. Probably, the high energy tails 
will break the scale-invariance at some stage and will 
close the system of moments by self-limitation of the tail 
of the distribution function. The more important require- 
ment here is that /(c p ,k) is a probability distribution 
function and must therefore be positive and real. This 
property renders 



/3/x/k < 1, 



< e p < oo 



(64) 



a condition that has been shown to be satisfied for ideal 
gases and in the absence of external potential fields. Oth- 
erwise one requires that 



> p. - k/P > 0, 



(65) 



where jl is the total chemical potential including the ex- 
ternal potential field. In the second case, a certain range 
of low particle energies (or low particle velocities) is not 
covered by the distribution. This implies that for suf- 
ficiently large positive total chemical potentials the mo- 
ment integrals may become non-analytic. Such cases pro- 
vide enormous difficulties for physical interpretation. In 
the following we will exclude this case from discussion. 

In terms of the velocity v of the particles the isotropic 
velocity distribution function reads 



/(v,«)=(l-^ + ^ 
K 2k 



-(k+1) 



(66) 



Since 1 — (3[i/k > 0, there exists a range of velocities for 
which 



mj3v 2 c /2K < 1 - Ph/k. 



(67) 



In this range of velocities the velocity distribution func- 
tion 



/(v,k)~(1-/Vk) 



-(re+l) _ 



const, v < v c (68) 



is practically flat, a property that has been frequently 
observed but was barely understood so far. The present 
theory provides a simple straightforward interpretation 
of this flatness problem. 

On the other hand, for high velocities v ^> v c the dis- 
tribution function becomes power law 



(69) 
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Any fit at high velocities with exponent a will provide 
a measurement of k — a/2 — 1. The break point of the 
distribution is at velocity 



v b » [2k(1 -/?/Vk)/to/?] 1/2 , 



(70) 



allowing for the determination of /3 = l/k B T, i.e. the 
determination of the temperature. 

A general important property of the (isotropic) veloc- 
ity distribution is that the value 4-7rp 2 /(p, k) is the prob- 
ability of finding any particles in the intervall of momen- 
tum [p,p + dp]. The maximum of this expression thus 
defined the most probable velocity v (momentum) of the 
particles. This velocity is given by 



2k B T 



1/2 r 



i qK T(k - 1/2) 



1/(2k-1) 



n T(k+1) 
On the other hand, the rms speed is defined as 



1 

N 



AttV 



d 3 pv 2 f(v,n) 



1/2 



The value obtained is 
3k B T^ 1/2 



K- 3/2 



i qR T(k - 1/2) 

n r(« + 1) 



(73) 



It is not surprising that this value differs from v as this 
is also true for the Maxwell-Boltzmann distribution of 
velocities. The ratio of both speeds is 



2 {k- 3/2 



(74) 



The most probable speed is smaller than the rms veloc- 
ity. The value of these expressions is mainly that mea- 
surement of both velocities provides an independent pos- 
sibility to determine the value of k and, subsequently, the 
value of the thermodynamic temperature T. 



B. Energy Distribution 

We briefly turn to a discussion of the energy distribution 
as this in practice is sometimes more important than the 
velocity distribution itself. Similar to the Boltzmann case 
the energy distribution is defined as 



(75) 



At small particle energies e it increases as / oc e 1 / 2 , 
Reaches maximum at 



f3e 1 — Ph/k 
~~K ~ 2k+1 ' 



(76) 



and, at high energies, becomes power-law and drops as 
e K+1 / 2 . Again, from its asymptotic behaviour one deter- 
mines k, while the value of the most probable energy is 
a measure of (3/1. One also notes that the most proba- 
ble energy increases with temperature T. The maximum 
value of the distribution itself is 



/m(e,«) 



(2 K m 3 ) 1 / 2 / 2k + 1 

[2( k + i)] k+1 yi-pn/K 



K+l/2 



(77) 



(71) 



(72) 



1/(2k-1) 



This value is always smaller than y/m 3 //3. Hence this 
value increases at about yT, Its determination allows 
for the measurement of the thermodynamic temperature 
of the gas and in combination with the above formula 
also for the measurement of the chemical potential. 

C. Gases in External Potentials 

An important frequently realised case is that a gas is 
brought into an external potential field <j>. Then the dis- 
tribution becomes 



/(e p -<£,«) = 



1 - 



Pi* . P(t P 



-(«+!) 



(78) 



The external potential simply changes the value of the 
chemical potential 



n + 4> 



(79) 



as in an ordinary gas as well. Depending on the sign of the 
external potential this may have enormous consequences 
on the analyticity of the distribution function, as has 
been mentioned in the introduction to this section. As 
long as this analyticity is guarantied, further integration 
of the distribution function provides no difficulties. 

As for a first and simple application we calculate the 
density of an electron plasma immersed into an external 
electric field. In this case the average external density in 
the absence of the field is no, and cj> — * — e<p in Equation 
(fr8|) is understood as the potential energy of the electron 
in the electric potential field <j). Carrying out the integra- 
tion of the zeroth order moment of [78] one finds for the 
electron density 



n{4>) 
no 



e<j) k - 1/2 
k B T k 



n Q r(«+i) 
n qK r(« - 1/2) 



1/(k-1/2) 



(80) 



This expression replaces the well-known commonly used 
Boltzmann formula 



Ub(4>) = n exp(e4>/k B T). 



(81) 



Clearly, this expression is the limiting form of Equation 
(p0|) for k — ► oo. An exactly equivalent formula may be 
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found, e.g., for gases embedded in an external gravita- 
tional field of acceleration — g (the so-called barometric 
law). In order to obtain the corresponding expression one 
replaces e<f>/k&T — > —mgH/k^T, where H is the height 
above level H = 0, where n(H = 0) = uq. The interest- 
ing point about both these density distributions is that, 
in a K-gas, the density reacts algebraically to the pres- 
ence of the external potential. For the barometric case 
this implies that it decays less than exponentially. 

A simple application of the density formula Equation 
( pp| ) is to the calculation of the screening distance of 
an ion in a quasineutral K-plasma, the so-called Debye 
length. This length may be found, expanding the right- 
hand side of Equation (80) for small potentials and us- 
ing Poisson's law. For illustration we show only the one- 
dimensional case. Then, from 



d 2 ^ 
dx 2 



e An 
So 



(82) 



(with G the vacuum dielectric constant) we find that 
the new Debye length Ad,k can be expressed through the 
Boltzmann-Debye length Adb as 



\ 2 - \ 2 

a D,k ~ A DB" 



3/2 



»,„r(K-i/2) 



n T(k + 1) 



1/0-1/2) 



(83) 



This value always exceeds the Boltzmann-Debye length. 
Intuitively this effect is clear and very satisfactory is as 
far as one may easily convince oneself that the presence 
of an excess of electrons in the tail of the distribution 
function over the Maxwell-Boltzmann distribution will 
worsen the conditions for screening. The faster electrons 
are less affected by the potential of the test ion and, 
hence, this potential will reach farther out into space. 
Actually, at the absolute limiting value K m i n = 3/2 the 
ratio of the two Debye-lengths becomes practically infi- 
nite, 



lim 



k^3/2 y a 



Ad,k 



>DB 



4 1 / 3 37rfc B T 
n^h 2 



3/4 



(84) 



The effect of such large screening lengths on the be- 
haviour of K-plasmas is of considerable interest. It may 
possibly support electric fields in plasmas stronger than 
ordinarily expected to exist. Also, because the level 
of electron thermal fluctuations in Langmuir waves is 
proportional to 1/noX 3 ^ K , Langmuir fluctuations in k- 
plasmas should be strongly reduced. It may be speculated 
that such plasmas would be even less collisional and ex- 
hibit lower levels of spontaneous emission as well as lower 
levels of thermal radiation from the plasma. However, no 
firm conclusion can be drawn at this level of investigation 
as long as the microscopic processes acting in K-plasmas 
have not been clarified. 



D. Relativistic Distribution 

In this last subsection we briefly investigate the relativis- 
tic version of the Loorentzian distribution function. In 
relativistic gases, e p = rwy^ic 2 , with 7 rc i the relativistic 
Gamma-factor. The phase space volume element trans- 
forms as 

Ai:d 3 xp 2 dp = 4tt V m 3 c 3 {^ cl - l) 1/2 d7 re i, (85) 

where dV — dVb/7 ro i,P = m-Troi = m c(7 2 c i — I) 1 / 2 . This 
leads to the following expression for the thermodynamic 
potential 



47r V m 3 c 3 



d7rd(7^-l) 1/2 



[1 - /3{i/k+ ((3mc 2 1 re)7 re i] 



(86) 



With the help of this expression we can calculate the 
particle number 



N = 



4ir Vbm 3 c 3 



d7rci(7 r 2 cl -l) 1/2 



[l-^/K+OSmc 2 /^^ 1 ' 



(87) 



From this formula we identify the relativistic distribution 
function as 



47rVbmV 2 1/2 

/rcl(7rol,K) = (7rcl ~ 1) ' 



0fi (3mc 2 

1 1 7rcl 

K K 



K > 2. 



The restriction on k results from the requirement that 
the relativistic energy as a moment of the distribution 
function must exist. For small but still relativistic ener- 
gies, i.e., for 

1 < 7rd < (K//3mc 2 )(l - Ph/k) (89) 

the relativistic distribution function varies according to 

/rcl(7rcl, K) - \/2( 7 rcl - 1) 1/2 . (90) 

This increase with energy reproduces the ^Je^ low en- 
ergy increase of the non-relativistic energy distribution 
function found above. On the other hand, in the ultra- 
relativistic domain 7 ro i ^> 1, 



/rel(7rel,«0 °C 7 re f> 



(91) 



implying a flatter decay than exhibited by the non- 
relativistic energy distribution, the latter evolving as 
f(e, k) cx e _ ( K + 1 / 2 ). The maximum of this distribution 
is at energy 



7rel 



1 — /3/i/k 
2(3mc 2 



1 



1 



k + 1 4/3 2 m 2 c 4 x 1/2 ' 
k 1 — /3/x/k 



(92) 



This is the most probable relativistic energy in a K-gas. 
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VII. DISCRETE ENERGY LEVELS: GIBBSIAN 
DISTRIBUTION 

The theory developed in 0] and explicated in the former 
sections of this paper assumes that the energy levels in 
the K-gas are continuously distributed. If, on the other 
hand, the energy levels are discrete, the formalism must 
be replaced by a discrete formalism. This can be achieved 
by replacing the integrals in the most fundamental for- 
mulas by sums over the discrete states of the system. 
The most important replacement is to be done in the 
thermodynamic potential Q of Equation (^). When sum- 
ming over discrete states, the volume factor V disappears, 
and the restriction on momentum space volume elements 
contained in the differential d 3 p/h 3 becomes unnecessary. 
Thus, the discrete partition function becomes 



if 

/3 



(3E r 



(93) 



Here E r is the r-th discrete energy level, and the sum- 
mation is over all dynamically possible levels E r . In a 
K-system, this expression replaces (up to an classically 
unimportant factor containing the chemical potential) 
the classical (Gibbsian) partition function 



JV 



(94) 



which is the sum over Gibbs factors, as can be shown 
by forming the limit value for k — > oo. (Note that the 
above mentioned factor —/3~ l exp(— /3/i) that remains af- 
terwards drops out by simple normalisation in the clas- 
sical case, while it is important to be retained in the 
K-system.) The partition function ( |93| ) contains all the 
information about the equilibrium state of the K-system. 
In analogy to our derivation of the distribution function 
it also suggests that the most probable Gibbs distribution 
of discrete states can be written as 



w r 



1 Pn | 0E r 

K K 



(95) 



Here w r is a probability (or else the most probable oc- 
cupation number of state r). Accordingly, the entropy of 
the discrete system becomes 



A" 



St = — ^b k 7 w r (l 



,-!/(«+!) 



), 



(96) 



an expression that is of similar though not identical form 
as the one given in ||. All thermodynamic relations ap- 
ply to Qt- It is important to note that there is no simple 
probabilistic way of distributing N particles on arbitrary 
sates that would reproduce the above Lorentz-Gibbsian 
probability distribution. The occupation number w r is 



not a simple most probable probability distribution in 
the ordinary meaning of the word. It cannot be achieved 
at by throwing the dice. Though it is obvious that it 
maximises St thereby leading to the state of maximum 
disorder, this does not mean that this state is achieved 
in a simple probabilistic way. There are non-obvious un- 
derlying processes that speed up this distribution. These 
processes will have to be investigated in future. 



VIII. CONCLUSIONS 

In the present paper we extended our initial theory of the 
statistical mechanics of the K-state to develop the ther- 
modynamics of such a state. We found that it is indeed 
possible and in asymptotic agreement with classical sta- 
tistical mechanics to construct the full thermodynamics 
of the K-state. It was possible to define the thermody- 
namic potential Qt that contains all the physics of the 
equilibrium state. One of the most important findings 
was that the definition of temperature can be retained in 
its full validity also in the K-state of a system. Another 
not less important conclusion was that the definition of 
observables in the K-state is exactly the same as in clas- 
sical statistical mechanics. This fact is most satisfactory 
because it does not violate physical intuition. Moreover, 
the entropy can only grow in states of the kind corre- 
sponding to K-states. This is satisfactory as well: addi- 
tion of disorder does not diminish disorder, it can only 
increase it. There is no process in closed systems, inde- 
pendent on their internal dynamics that could minimise 
entropy. To this end one must do external work on the 
system. The claim that instability would generate order 
in a closed system and reduce entropy is misleading. En- 
tropy can only be reduced locally on the expense of en- 
tropy increase in other locations in the system. Instability 
can do nothing else but redistribute free energy until it 
makes it available to other dissipative processes. 

A large number of questions has been left open for 
future investigation. The most important is not that 
concerning the reality of processes like the once dis- 
cussed here. Observation of Levy flights and increasing 
frequency of measurement of so-called '^-distributions' 
(after those we have modelled our distribution function 
by the use of the control parameter n) in large collision- 
less systems like the near-Earth space environment pro- 
vide sufficient justification for an expedition into basic 
physics. But the question remains of what is the phys- 
ical relevance of the index k! In an interesting paper 
p5|, Hasegawa and co-workers found a similar distribu- 
tion function solving a Fokker-Planck equation under the 
assumption of a very particular weakly turbulent inter- 
action in a plasma. There, k turned out to depend on the 
dispersion of plasma waves. 

A solution of this kind may remain a particular case 
of weak turbulence. It rather seems that the k state of a 
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system refers to an intermediate strongly turbulent state 
when the system has settled into quasi-equilibrium for 
some time before binary collisions set on to destroy the 
state and to dissipate the energy stored in the turbu- 
lent motion, as has been argued in |Ej . Equilibria usually 
imply the dissipative destruction of the free energy. In n- 
cquilibrium it seems that this is not necessarily the case. 
But the observation that k appears only in this inter- 
mediate state, suggests that there must be discontinuous 
transitions between the initial and the K-state and pos- 
sibly also to the Boltzmann regime. 

Figure |l| schematically illustrates this behaviour of k 
for the same scenario as in |Q. In this figure we plot- 
ted the value of k~ 1 . This value must be zero in the 
respective linear and in the final Boltzmann regimes, be- 
cause these regimes are classical (k — * oo). In the tur- 
bulent quasi-equilibrium state n~ l suddenly assumes a 
value between zero and 1/K m in (the figure shows the 
marginal case k = const = K m in)- The sudden increase 
from kT 1 = to n^ 1 ^ is understood as a critical phase 
transition. A possible non-critical transition is shown as 
the dotted line in the non-linear violent relaxation regime 
where our theory does not apply. 

Physically spoken, there is no need for k being con- 
stant troughout the entire quasi-equilibrium. If k changes 
sufficiently slowly compared to the length of this quasi- 
equilibrium the equilibrium theory can still be applied. 
There is reason to believe that such a situation is closer 
to reality than k = const. Dissipative non-collisional ir- 
reversible processes will necessarily accompany the tur- 
bulent state. Because the heating causes T to increase 
with time, k(T) will itself become a function of time. 
Two such cases have been included in Figure |. In the 
first one (thin solid line), the final critical phase transi- 
tion does entirely disappear. k(T) tends to infinity right 
at the end of the intermediate stationary state. In the 
second case (dotted line), k remains finite towards the 
end of the quasi-equilibrium, and another critical phase 
transition is needed when suddenly binary collisions take 
over. Which one of these cases will be realised, depends 
on the properties of the system during its evolution. We 
will have to develop the non-equilibrium thermodynam- 
ics for K-systems before these questions can be answered. 
This goal in addition requires to solve the problem of 
divergencies appearing in the higher moments of the dis- 
tribution function. 
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FIG. 1. Schematic evolution of 1/k for a different cases. 
The solid line refers to the model used in the present paper. 
k is zero At times t < t n i and t > t c = l/nu c , respectively. 
There are two critical transitions in this case where k jumps 
from zero to a finite value k > K m in at these transitions (the 
figure shows the marginal case k — /t m i n )- The dotted line 
shows a possible Non critical transition where k ^ already 
in the interval ti < t < i„; . Such a transition may avoid criti- 
cality and would be describable by ordinary strong turbulence 
theory (nonlinear wave- wave and wave-particle interactions). 
In addition, k can depend on time (or otherwise temperature 
if accounting for transitional heating in turbulence. Then ei- 
ther the cases of the dotted (retaining a late critical phase 
transition) or thin solid lines (no final critical transition) ap- 
ply, respectively. 
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